Isospin Dependence of 1 5'o Proton and Neutron Superfluidity in Asymmetric Nuclear 

Matter 
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We investigate the 1 So neutron and proton superfluidity in isospin asymmetric nuclear matter. We 
have concentrated on the isospin dependence of the pairing gaps and the effect of a microscopic three- 
body force. It is found that as the isospin asymmetry goes higher, the neutron So superfluid phase 
shrinks gradually to a smaller density domain, while the proton one extends rapidly to a much wider 
density domain. The three-body force turns out to weaken the neutron 1 So superfluidity slightly, 
but it suppresses strongly the proton So superfluidity at high densities in nuclear matter with large 
isospin asymmetry. 
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Superfluidity plays an important role in understand- 
ing a number of astrophysical phenomena in neutron 
stars Q S S S H B E S In the inner crust of 
a neutron star where the total baryon density p is low, 
neutrons are thought to be superfluid in the 1 Sq chan- 
nel. In the nuclear core part (i.e, the outer core part), 
protons may form a superfluid in the partial wave 
states, and neutron superfluidity is expected in the P2- 
3 F2 partial wave channel. It is generally expected that 
the cooling processes via neutrino emission S 1 the 
properties of rotating dynamics, the post-glitch timing 
observations 0, Q, the possible vertex pinning || of a 
neutron star are rather sensitive to the presence of neu- 
tron and proton superfluid phases as well as to their pair- 
ing strength. Recently it has been shown ; 9| that a weak 
proton 1 jS'o superfluidity is in consistent with the cool- 
ing data if the effect of the possible presence of accreted 
envelopes is taken into account in the cooling scenario. 

Nucleon superfluidity in symmetric nuclear matter and 
pure neutron matter has been investigated extensively 



All these investigations have predicted the occurrence of 
the So nucleon superfluid phase in the low density re- 
gion although the obtained strengthes of the superfluid- 
ity are remarkably sensitive to the different approaches 
and the different approximations adopted. There also 
exists few investigations on the nucleon superfluidity in 
/3-stable neutron star matter [25| . It is shown that due to 
the small proton fraction in /3-stable matter, the proton 
superfluid phase in the 1 So channel may extend to much 
higher baryon densities with a much smaller maximum 
of the pairing gap as compared to the case of symmetric 
nuclear matter. Therefore it is of interest to investigate 



the variation of nucleon superfluidity vs. isospin asym- 
metry in asymmetric nuclear matter, which is expected 
to be helpful for understanding the properties of nucleon 
superfluid phases in neutron stars. 

The aim of this work is devoted to the isospin depen- 
dence of the neutron and proton superfluid phases 
in asymmetric nuclear matter, and to investigating the 
influence of three-body forces which turn out to be cru- 
cial for reproducing the empirical saturation properties 
of nuclear matter in a non-relativistic microscopic ap- 
proach [2rll27iE^ |. 

For such a purpose, we shall not go beyond the BCS 
framework. In this case, the pairing gap which character- 
izes the superfluidity in a homogeneous Fermi syst em is 
determined by the standard BCS gap equation |22j, i.e., 



A £ = -X>(M')^7A £ , 



(1) 



'Corresponding address: Institute of Modern Physics, Chinese 
Academy of Science, P.O. Box 31, Lanzhou 730000, P.R. China. 
Tel: 0086-931-4969318; E-mail: zuowei@impcas.ac.cn 



where v(k, k') is the realistic NN interaction in momen- 
tum space, Ej: = ^/(eg — £f) 2 + A|, eg and €f being the 

s.p. energy and its value at the Fermi surface, respec- 
tively. In the BCS gap equation, the most important in- 
gredients are the realistic NN interaction v(k, k') and the 
neutron and proton s.p. energies eg in asymmetric nu- 
clear matter. For the NN interaction, we adopt the Ar- 
gonne V18 (AVis) two-body interaction [3(| plus a micro- 
scopic three-body force (TBF). The TBF adopted in the 
present calculation was originally proposed in Ref. j2|| 
based on the meson-exchange current approach. The pa- 
rameters of the TBF, i.e., the coupling constants and the 
form factors, were determined |27j from the one-boson- 
exchange potential model to meet the self-consistent re- 
quirement with the adopted AVis two-body force. A 
more detailed description of the TBF model and the re- 
lated approximations can be found in Ref. |2pj |. 

The proton and neutron s.p. energies in asymmetric 
nuclear matter are calculated by using the BHF approach 
for isospin asymmetric nuclear matter |3lj |. The starting 
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point of the BHF approach is the Brueckner G-matrix, 
which satisfies the Bethe-Goldstone (BG) equation |3^ |: 



kik2 



\kik 2 )Q(k 1 ,k 2 )(k 1 k 2 \ 
uj - e(fci) - e(fc 2 ) + it] 



G(p,ft;u), 
(2) 

where ki = (fci,<7j,Tj) denotes the momentum, the z- 
components of spin and isospin of a nuclcon, respectively. 
u) is the starting energy and Q(fci, k 2 ) is the Pauli oper- 
ator. The isospin asymmetry parameter is defined as 
ft = (pn — Pp)l Pi being p n , p p , and p the neutron, proton 
and total nucleon number densities, respectively. The s.p. 
energy is given by e(fc) = h 2 k 2 /(2m) + U(k). In solving 
the BG equation, we adopt the continuous choice [32( for 
the s.p. potential U(k) since it has been proved to pro- 
vide a much faster convergency than the gap choice [23 • 
Under the continuous choice, the s.p. potential describes 
physically the nuclear mean field felt by a nucleon in nu- 
clear medium |3j| and is calculated from the real part of 
the on-shell G-matrix. 

The neutron and proton superfluidity in nuclear mat- 
ter can be described by their pairing energy gaps at their 
respective Fermi surfaces. To solve the gap equation, we 
follow the scheme given in Ref. |To| . The proton and neu- 
tron s.p. energies in the gap equation are calculated from 
the BHF approach. In our calculations, the TBF contri- 
bution has been included by reducing the TBF to an 
equivalently effective two-body interaction according to 
the standard scheme as described in Ref. [2(| . A detailed 
description and justification of the method are discussed 
in Refs. [H, E3 

In Fig. H is shown the neutron energy gap in the 1 S'o 
partial wave channel = A(kp) as a function of the 
total baryon density p. The curves along the direction of 
the arrow from the bottom correspond to (3 = 0.2, 0.4, 
0.6, and 0.8, respectively. We see that the neutron 1 Sq 
superfiuid phase exists only at low densities (p < 0.13 
fm -3 ) and the peaks of the pairing gaps are located 
around p = 0.02 fm~ 3 , which is compatible with the pre- 
vious predictions for pure neutron matter and symmet- 
ric nuclear matter [101 l2-l| . As the isospin asymmetry ft 
increases, the neutron fraction increases for a given to- 
tal density and as a consequence the density domain for 
the neutron So superfluidity shrinks gradually. In the 
case without including the TBF, the density region for 
the presence of the neutron superfluidity reduces from 
p < 0.13 fm~ 3 down to p < 0.1 fm~ 3 as the isospin 
asymmetry rises from /3 = 0.2to/3 = 0.8. In Ref.[24j, 
we calculated the 1 Sq pairing gaps in symmetric nuclear 
matter (ft = 0) and pure neutron matter (ft = 1). We 
found that the neutron 1 Sq superfiuid phase is expected 
to appear in a density range of ps < 0.135 fm~ 3 for 
symmetric nuclear matter and pb < 0.09 fm~ 3 for pure 
neutron matter, in consistent with our present results. 
The above obtained isospin dependence of the pairing 
gap is readily understood as follows. In asymmetric nu- 



clear matter, the strength of the neutron pairing gap is 
related directly with the neutron number density. As 
the asymmetry rises, the neutron excess and the neu- 
tron number density for a given total nucleon density 
increase. As a consequences, the density domain for the 
existence of the neutron superfiuid phase shrinks gradu- 
ally going from symmetric matter to pure neutron mat- 
ter. We notice that in both cases with and without in- 
cluding the TBF, the peak values of the gaps decrease 
and their location shift slightly to lower densities as the 
matter becomes more isospin asymmetric. This is mainly 
attributed to the isospin dependence of the neutron s.p. 
energy spectrum in asymmetric nuclear matter. As we 
know, the neutron s.p. potential in asymmetric matter 
becomes less deeper and the neutron effective mass gets 
larger for a larger isospin asymmetry according to the 
BHF and DBHF calculations |3l[, |35| . We may also see 
from Fig. ^ that within the region of the neutron su- 
perfiuid phase, as the isospin asymmetry increases, the 
pairing gap becomes smaller at low densities below and 
in the vicinity of the peak position, while it gets larger at 
relatively high densities. The above obtained isospin be- 
havior stems from the competition between two different 
mechanisms. On the one hand, at a fixed total density 
on the right side of the peak position, a higher asymme- 
try corresponds to a larger value of neutron density and 
a weaker neutron superfluidity. On the other hand, the 
isospin dependence of the neutron s.p. spectrum intends 
to enhance the neutron pairing gap as the isospin asym- 
metry increases. At low densities around and below the 
peak position, the later effect is dominated, while as the 
density increases, the former mechanism becomes more 
and more effective. Comparing the solid curves with the 
corresponding dashed ones, we may see that the TBF af- 
fects mainly the pairing gap at high densities. Inclusion 
of the TBF weakens the neutron superfluidity at high 
densities and makes the predicted density range for the 
existence of neutron superfiuid phase smaller. 

In Fig. [21 is reported the proton 1 Sq pairing gap in 
asymmetric nuclear matter for four different asymmetries 
ft = 0.2, 0.4, 0.6, and 0.8, respectively. The solid curves 
arc obtained by including the TBF, while the dashed ones 
without including the TBF. The isospin dependence of 
the proton 1 Sq pairing gaps turns out to be completely 
different from that of the corresponding neutron ones in 
both cases with and without the TBF. It is seen that 
as the isospin asymmetry increases, the density domains 
for the existence of the proton superfiuid phases enlarge 
rapidly, the peaks of the pairing gaps become lower ap- 
preciably and shift to higher densities gradually. As the 
asymmetry rises from ft — 0.2 to ft = 0.8, the density 
domain extends from p < 0.17 fm~ 3 to p < 0.435 fm -3 
in the case of not including the TBF and from p < 0.13 
fm~ 3 to p < 0.25 fm~ 3 in the case of including the TBF, 
the peak values lower from about 2.0 MeV to about 1.4 
MeV in both cases. The difference between the proton 
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and neutron superfluid phases is especially pronounced 
at high isospin asymmetries. As compared to the neu- 
tron 1 Sq superfluidity, the proton 1 S superfluid phase in 
highly asymmetric matter extends to much higher den- 
sities but with a smaller peak value of the pairing gap. 
For example, in the case without including the TBF, at 
(3 = 0.8, the proton 1 Sq superfluidity is predicted to exist 
in the density region of p < 0.435 fm~ 3 which is much 
wider than that of p < 0.1 fm~ 3 for the neutron 1 So su- 
perfluidity The above behavior of the proton superflu- 
idity vs. isospin asymmetry can be readily explained in 
terms of the isospin dependence of the proton density and 
the proton s.p. spectrum in asymmetric nuclear matter. 
First, at a fixed total density, a higher asymmetry cor- 
responds to a smaller proton concentration and a lower 
proton density. Such an isospin dependence of the proton 
density is directly responsible for the widening of the den- 
sity domain as a function of asymmetry. Second, the pro- 
ton s.p. potential becomes deeper going from symmetric 
matter to pure neutron matter [3ll| . and this results in 
the lowering of the proton pairing peak vs. asymmetry. 
At relatively low densities below and around the peak 
position, the variation of the proton s.p. spectrum as a 
function of asymmetry plays an major role in determin- 
ing the isospin variation of the proton pairing gap and 
thus the proton superfluidity becomes weaker at a higher 
asymmetry. While as the density increases, the effect due 
to the decreasing of the proton fraction vs. asymmetry 
gets stronger. At high enough densities, it becomes pre- 
dominant and leads to a turnover of the isospin behavior 
of the proton superfluidity (i.e., a stronger superfluidity 
at a higher asymmetry). The turnover is clearly seen in 
Fig. 121 and stems from the competition between the two 
above-mentioned isospin effects. 

We notice that the TBF effect is negligibly small at 
low densities below and around the peak position. How- 
ever, it gets stronger rapidly as the density goes up. The 
TBF turns out to induce a significant reduction of the 
proton pairing gaps in the high-density superfluidity do- 
main. Consequently, it leads to a remarkable shrinking 
of the density domain for the existence of the proton 
superfluid phase. We notice that the above predicted 
TBF suppression of the 1 Sq proton superfluidity is par- 
ticularly pronounced for highly asymmetric matter. For 
example, at j3 = 0.8, the density domain of the pro- 
ton superfluid phase is reduced by about 50%, i.e., from 
p < 0.435 fm -3 down to p < 0.25 fm -3 , by inclusion 
of the TBF. In nuclear matter, proton pairs are embed- 
ded in the medium of neutrons and protons, both the 
surrounded protons and neutrons contribute to the TBF 
renormalization of the proton-proton interaction, there- 
fore the TBF effect on the proton pairing correlations is 
determined by the total nucleon number density instead 
of the proton number density. Accordingly, in spite of 
the small proton fractions and the low proton densities at 
high asymmetries, the TBF modifies strongly the proton- 



proton pairing interactions at high-density superfluidity 
domain and weakens considerably the corresponding pro- 
ton pairing gaps. One can verify readily from Fig. [3 
that the TBF suppression of the proton superfluidity is 
mainly in the high-density region and the reduction of the 
gap increases rapidly as increasing the total density. We 
also notice from Fig.[3]that the TBF effect gets stronger 
as the asymmetry increases since the proton superfluid 
phase extends to larger densities for asymmetric mat- 
ter at higher isospin asymmetries. Inclusion of the TBF 
weakens considerably the isospin dependence of the pre- 
dicted proton superfluidity as compared to the results 
obtained by adopting purely the AVi$ two-body force. 

In summary, we have calculated the neutron and pro- 
ton 1 So pairing gaps in asymmetric nuclear matter based 
on the BHF approach and the BCS theory. We have espe- 
cially investigated the isospin dependence of the pairing 
gaps in the 1 S'o channel and the influence of the TBF. It 
is shown that the isospin dependence of the proton 1 Sq 
superfluidity in asymmetric nuclear matter is completely 
different from that of the neutron one. The neutron 1 Sq 
superfluid phase exists only in low density region for all 
isospin asymmetries. As the matter goes from symmetric 
nuclear matter to pure neutron matter, the peak value 
of the neuron 1 Sq pairing gap becomes larger and the 
density domain for the existence of neutron superfluidity 
shrinks gradually. While the density domain for the pro- 
ton superfluid phase enlarges rapidly as the isospin asym- 
metry rises and it may extend to very high densities for 
highly asymmetric nuclear matter. The peak value of the 
proton pairing gap turns out to be larger in asymmetric 
nuclear matter at a higher isospin asymmetry. 

The TBF affects only weakly the neutron 1 So super- 
fluid phase in asymmetric nuclear matter, i.e., it reduces 
slightly the pairing gap, due to the low-density region for 
this kind of superfluidity. However it suppresses strongly 
the proton superfluidity in the 1 Sq channel at high densi- 
ties, especially at high asymmetries. The density domain 
for the existence of the proton 1 So superfluid phase is re- 
duced by about 50% from p < 0.435 fm~ 3 to p < 0.25 
fm~ 3 by inclusion of the TBF. 
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Commission, and DFG, Germany. 



[1] D. J. Dean and M. Hjorth- Jensen, Rev. Mod. Phys. 75, 

607 (2003) and references therein. 
[2] U. Lombardo and H.-J. Schulze, Lecture Notes in 

Physics, vol. 578, Eds. D. Blaschke, N. K. Glendenning 

and A. Sedrakian (Springer- Verlag, Berlin, 2001), p. 30. 
[3] D. Page, Astrophys. J. 428, 250 (1994). 
[4] H. Heiselberg and M. Hjorth- Jensen, Phys. Rep. 328, 237 

(2000). 



4 



[5] B. Link, Phys. Rev. Lett. 91, 101101(2003). 

[6] S. L. Shapiro and S. A. Teukosky, Black Holes, White 

Dwarfs and Neutron Stars (Wiley, New York, 1983). 
[7] J. A. Sauls, Timing Neutron Stars, eds. H. Ogelman and 

E. P. J. van den Henvel, (Dordrecht, Kluwer, 1989) p.457. 
[8] D. Pines and M. A. Alpar, Nature, 316, 27 (1985). 
[9] M. E. Gusakov, A. D. Kaminker, D. G. Yakovlev, and O. 

Y. Gnedin, Mon. Not. R. Astron. Soc. 363, 555 (2005); 

A. D. Kaminker, M. E. Gusakov, D. G. Yakovlev, and O. 

Y. Gnedin, Mon. Not. R. Astron. Soc. 365, 1300 (2006). 
[10] M. Baldo, J. Cugnon, A. Lejeune and U. Lombardo, 

Nucl. Phys. A515, 409 (1990). 
[11] 0. Elgaroy, L. Engvik, M. Hjorth- Jensen and E. Osnes, 

Nucl. Phys. A604, 466 (1996). 
[12] T. Duguet, Phys. Rev. C69, 054317 (2004). 
[13] A. Fabrocini, S. Fantoni, A. Y. Illarionov, and K. E. 

Schmidt, Phys. Rev. Lett. 95, 192501 (2005). 
[14] J. W. Clark, C. G. Kallman, C.-H. Yang and D. A. 

Chakkalakal, Phys. Lett. B61, 331 (1976); J. M. C. Chen, 

J. W. Clark, R. D. Dave and V. V. Khodel, Nucl. Phys. 

A555, 59 (1993). 
[15] T. L. Anisworth, J. Wambach and D. Pines, Phys. Lett. 

B222, 173 (1989); J. Wambach, T. L. Anisworth and 

D. Pines, Nucl. Phys. A555, 128 (1993). 
[16] C. J. Pethick and D. G. Rawenhall, Ann. N. Y. Acad. Sci. 

647, 503 (1991). 
[17] H.-J. Schulze, J. Cugnon, A. Lejeune, M. Baldo and 

U. Lombardo, Phys. Lett. B375, 1 (1996). 
[18] M. Serra, A. Rummel and P. Ring, Phys. Rev. C65, 

014304 (2001). 

[19] A. Rabhi, R. Bennaceur, G. Chanlray and P. Schuck, 

Phys. Rev. C66, 064315 (2002). 
[20] A. Schwenk, B. Friman and G. E. Brown, Nucl. Phys. 

A713, 191 (2003). 
[21] H. Muther and W. H. Dickhoff, Phys. Rev. C72, 054313 

(2005). 

[22] M. Baldo and A. Grasso, Phys. Lett. B485, 155 (2000); 

U. Lombardo, P. Schuck and W. Zuo, Phys. Rev. C64, 

021301(R), (2001). 
[23] P. Bozek, Phys. Rev. C62, 054316 (2000); Phys. Lett. 

B551, 93 (2003). 
[24] W. Zuo, U. Lombardo, H. J. Schulze, and C. W. Shen, 

Phys. Rev. C66, 037303, (2002). 
[25] M. Baldo, J. Cugnon, A. Lejeune and U. Lombardo, 



Nucl. Phys. A536, 349 (1992); 0. Elgaroy, L. Engvik, 

M. Hjorth- Jensen and E. Osnes, Phys. Rev. Lett. 77, 

1428 (1996); W. Zuo, Z. H. Li, G. C. Lu et al, Phys. 

Lett. B595, 44 (2004); T. Takatsuka and R. Tamagaki, 

Nucl. Phys. A738, 387 (2004). 
[26] P. Grange, A. Lejeune, M. Martzolff and J.-F. Mathiot, 

Phys. Rev. C40, 1040 (1989). 
[27] W. Zuo, A. Lejeune, U. Lombardo and J.-F. Mathiot, 

Nucl. Phys. A706, 418 (2002); Eur. Phys. J. A14, 469 

(2002). 

[28] R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989); 
M. Baldo, I. Bombaci, and G. F.Burgio, Astronomy 
and Astrophysics 328, 274 (1997); B. S. Pudliner, 
V. R.Pandharipande, J. Carlson, and R. B. Wiringa, 
Phys. Rev. Lett. 74, 4396 (1995). 

[29] P. Ring and P. Schuck, The Nuclear Many Body Problem 
(Springer- Verlag, New York, 1980). 

[30] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, 
Phys. Rev. C51, 38 (1995). 

[31] W. Zuo, I. Bombaci and U. Lombardo, Phys. Rev. C60, 
024605, (1999); W. Zuo, L. G. Cao, B. A. Li, U. Lom- 
bardo and C. W. Shen, Phys. Rev. C72, 014005 (2005). 

[32] J. P. Jeukenne, A. Lejeune and C. Mahaux, Phys. Rep. 
25C, 83 (1976). 

[33] H. Q. Song, M. Baldo, G. Giansiracusa, and U. Lom- 
bardo, Phys. Rev. Lett. 81, 1584 (1998). 

[34] A. Lejeune, C. Mahaux, Nucl. Phys. A 295, 189 (1978). 

[35] E. N. E. van Dalen, C. Fuchs and A. Faessler, Phys. Rev. 
Lett. 95, 022302 (2005); Phys. Rev. C72, 065803 (2005). 

Figure Captions: 



FIG. 1: Neutron 1 So energy gap as a lunction of density in 
asymmetric nuclear matter at various asymmetries. The solid 
curves are predicted by adopting purely the A Vis two-body 
interaction, while the dashed ones by using the AVi& plus the 
TBF. 



FIG. 2: The same as Fig.Qbut for proton 1 5'o pairing gap. 
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